Non‐intrusive reduced‐order modeling for multiphase porous media flows using Smolyak sparse grids

In this article, we describe a non‐intrusive reduction method for porous media multiphase flows using Smolyak sparse grids. This is the first attempt at applying such an non‐intrusive reduced‐order modelling (NIROM) based on Smolyak sparse grids to porous media multiphase flows. The advantage of this NIROM for porous media multiphase flows resides in that its non‐intrusiveness, which means it does not require modifications to the source code of full model. Another novelty is that it uses Smolyak sparse grids to construct a set of hypersurfaces representing the reduced‐porous media multiphase problem. This NIROM is implemented under the framework of an unstructured mesh control volume finite element multiphase model. Numerical examples show that the NIROM accuracy relative to the high‐fidelity model is maintained, whilst the computational cost is reduced by several orders of magnitude. Copyright © 2016 John Wiley & Sons, Ltd.     

The Impact of Tides on the Capillary Transition Zone

The capillary transition zone, also known as the capillary fringe, is a zone where water saturations decrease with height above the water table/oil–water contact as a result of capillary action. In some oil reservoirs, this zone may contain a significant proportion of the oil in place. In groundwater assessments, the capillary fringe can profoundly affect contaminant transport. In this study, we investigated the influence of a tidally induced, semi-diurnal, change in water table depth on the water saturation distribution in the capillary fringe/transition zone. The investigation used a mixture of laboratory experiments, in which the change in saturation with depth was monitored over a period of 90 days, and numerical simulation. We show that tidal changes in water table depth can significantly alter the vertical water saturation profile from what would be predicted using capillary–gravity equilibrium and the drainage or imbibition capillary pressure curves.     

The Effect of Viscosity Ratio and Peclet Number on Miscible Viscous Fingering in a Hele-Shaw Cell: A Combined Numerical and Experimental Study

Transport in Porous Media - The results from a series of well characterised, unstable, miscible displacement experiments in a Hele-Shaw cell with a quarter five-spot source-sink geometry are...     

Predicting Vertical Flow Barriers Using Tracer Diffusion in Partially Saturated,Layered Porous Media

Sudden changes in isotopic tracer concentration in pore waters have been interpreted as indicating barriers to vertical advective flow through porous rocks in the subsurface, e.g. step changes in \(^{87}\hbox {Sr}/^{86}\) Sr ratio are often used in the oil and gas industry as a signature of reservoir compartmentalisation. This study shows that this is not necessarily the case. It can take millions of years for such step changes to equilibrate by diffusion if there is no flow resulting from pressure or density gradients even in high permeability, high porosity rocks, particularly if the water saturation is low. Changes in tracer concentration gradients can be good indicators of changes in porosity (or water saturation) between layers. In contrast changes in sorption without a change in porosity are almost impossible to identify. The time taken for concentration gradients to equilibrate is affected by the layer properties but can be quickly estimated from the harmonic average of the effective diffusion coefficient for each layer and a simple analytical expression for a homogeneous system. This was achieved by performing a sensitivity analysis on different layer properties (porosity contrast, saturation contrast, sorption contrast, thickness ratio) using existing analytical solutions for diffusion in layered systems.     

Intrinsic equations for a generalized relaxed elastic line on an oriented surface in the galilean space

In this article, we derive the intrinsic equations for a generalized relaxed elastic line on an oriented surface in the Galilean 3-dimensional space G₃. These equations will give direct and more geometric approach to questions concerning about generalized relaxed elastic lines on an oriented surface in G₃.